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Abstract 

We consider a system of reaction-diffusion equations describing the 
reversible reaction of two species U, V forming a third species W and vice 
versa according to mass action law kinetics with arbitrary stochiometric 
coefficients (equal or larger than one). 

Firstly, we prove existence of global classical solutions via improved 
duality estimates under the assumption that one of the diffusion coeffi¬ 
cients of W or V is sufficiently close to the diffusion coefficient of W. 

Secondly, we derive an entropy entropy-dissipation estimate, that is a 
functional inequality, which applied to global solutions of these reaction- 
diffusion system proves exponential convergence to equilibrium with ex¬ 
plicit rates and constants. 
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1 Introduction 


The aim of this paper is to investigate the following class of systems of three 
nonlinear reaction-diffusion equations. For stoichiometric coefficients a, /3 ,7 > 
1 , we consider the system 


dtu — diAxU =— kw'^), (1) 

dtv — d 2 AxV =— kw'^), ( 2 ) 

dtw — dsAxW = 'y{£ u°‘v^ — k w'^), (3) 

where the constants i,k > 0 are positive reaction rates and di, ^ 2 , da > 0 denote 
positive diffnsion coefficients. Moreover, we assume that the concentrations u, v 
and w satisfy homogeneous Neumann conditions 

n{x)-yxU = Q, n{x)-WxV = Q, n{x)-VxW = 0 x G (4) 

on a bounded domain 17 with sufficiently smooth boundary 917 and the non¬ 
negative initial condition 

m( 0, a:) = uo(x) > 0, u(0, a:) = uo(a;) > 0, a(;(0, a;) = a(;o(a:) > 0 . (5) 

We remark that the reaction dynamics of System (l)-(3) is positivity pre¬ 
serving (see e.g. [15]) and, thus, that non-negative initial data (5) ensures 
non-negative solutions u, v, and w, i.e. 

(P) u{t,x) > 0, v{t,x) > 0, w{t,x) >0, Vt > 0, a: € 17. 

Moreover, the System (l)-(3) with homogeneous Neumann boundary con¬ 
ditions (4) satisfies the following two mass conservation laws: By taking the 
sums 7 X (1) -I- a X (3) and 7 x (2) -|- /3 x (3) and by integrating by parts with 
non-flux boundary conditions, we obtain the following conservation laws to hold 
for solutions of (1) - (5): 


/ {-fu{t,x) + aw{t,x)) dx = / ( 7 Uo(a:) -I- q:wo(x)) da: =: Ml > 0, 
Jq, Jq 


/ {'-fv{t,x) +/3w{t,x)) dx 
JQ 


{■jvo{x) + /3wo{x)) dx =: M 2 > 0, 


( 6 ) 


where we assume strictly positive masses Mi, M 2 > 0 . 

We remark that the System (l)-(3) can be further simplified by rescaling 
times, space and concentrations. In particular for all stoichiometric coefficients 
a -I- /3 ^ 7 , the rescaling 


(I) 


^1 , 
k ’ 


|17|' 






{u,v,w), 


allows (without loss of generality) to consider the System (1) - (3) with nor¬ 
malised reaction rates £ = fc = 1 on a normalised domain, i.e. we can replace in 
all three equations (1) - (3) 


{£u°‘v^ — k w^) —>■ (u“u^ — w'^), |f7| = 1. 
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In the special cases when a + (3 = 7 , it is still possible, for instance, to 
rescaling {£u°‘v^ — k w'^) —>■ — w'^) in the first two equations ( 1 ) and ( 2 ), 

but the third equation (3) can only be rescaled up to a common factor 
i.e. ^ 

"f{£u°‘v^ — k w'^) —?► 7 ^ ^ [u°‘v^ — vf^). 

For the rest of the article, we shall always consider the complete rescaled 
system 


dtu — diAxU = —a{u°‘v^ — w'^), 
dtv — d 2 AxV = —/3{u°‘v^ — w'*'), 
dtw — dsA^w = 'y{u°‘v^ — w^), 


(7) 

( 8 ) 
(9) 


subject to the initial date (5) and the homogeneous Neumann boundary con¬ 
ditions (4). In the spacial cases a -|- /3 = 7 , the additional factor in (9) 

leads to non-essential modifications of the arguments and calculations in this 
paper. 

The System (7) - (9) with the conservation laws ( 6 ) possesses a unique 
detailed balance equilibrium state, which are the unique positive constants 
{uoo,Voc,Woo) balancing the three reversible reactions, i.e. 

^00 ^00 ^00 

and satisfying the conservation laws ( 6 ), i.e. 

7U00 + awoo = Ml, ^Voo + / 3 woo = M2 

(recall |n| = 1). In fact, the uniqueness of the equilibrium {uoc,Voc,Woc) follows 
from the uniqueness of the solution Wao of the transcendental equation 


Ml _ Wooa\ f ^ _ Woo/3 ' _ 

7 1 J \ 1 7 


( 10 ) 


where the left-hand-side of ( 10 ) is strictly monotone decreasing in Woo while the 
right-hand-side of ( 10 ) is strictly monotone increasing in Woo- 

Finally, the System (7) - (9) features a Boltzmann-type entropy functional 
(or physically, a free energy functional) 

E{u,v,w){t) = [ (M(ln(u) — 1)-I-i;(ln(r;) — 1)-I-r(;(ln(u>) — 1)) dx, (11) 
Jn 

which dissipates according to the non-negative entropy dissipation functional 

O = -iE- 


Jn u Jn V 


+ / ( 

Jn 




In w 


,„(lf^),,>0. (12) 
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We remark, that thanks to the properties of the equilibrium {uoo,Voo,Woo), one 
can rewrite the relative entropy with respect to the equilibrium state in the 
following way: 


E{u,V,w){t) - E{u^,Voo,Woo) = / ( Ulnf — ) - (u - Moo) 

Jn \ 

+Mln( — ) - {v - Moo) + wln() - {w - Woo) ) dx, (13) 
VMoo/ ^Woo-' ) 

Clearly, the relative entropy dissipates according to the same entropy dissipation 
functional, i.e. 

^E{u,v,w){t) = ■^{E{u,v,w) - E{uoo,Voo,u)oo)){t) = -D{u,v,w){t). 

at at 

Remark 1.1. We emphasise that the entropy dissipation functional D(u,v,w) 
vanishes for all constant states (u,v,w)=(u,v,W), which balance the reactions, 
i.e. These states, however, form a two parameter family and 

one crucially requires the two conservation laws ( 6 ) as additional constraints 
in order to identify the equilibrium (uoo,Voo,Woo) os the unique state, which 
minimises/annihilates the entropy dissipation functional along the flow of mass- 
conserving solutions. 

The following Theorem is the first main result of this paper. The proof is 
based on applying (point-wise in time) a functional inequality between relative 
entropy and entropy-dissipation, see Proposition 3.1 below : 

Theorem 1.1. Let Lt be a connected open set of {N > 1) with sufficiently 
smooth boundary dQ such that Poincare’s inequality and the Logarithmic Sobolev 
inequality hold. Let c?i,c? 2 ,d 3 > 0 be three strictly positive diffusion coefficients 
and let a, fi,"f > 1. 

Consider non-negative initial data {uo,vo,wo) with finite mass and finite 
entropy E(uo,vo,wo) and assume that {u,v,w) is a non-negative solution (be it 
either classical, weak or possibly even renormalised, see Remark 1.2) of System 
(7) - (9) on an time interval 0 < t < T for T < -|-oo, which is assumed 

i) to satisfy the conservation laws ( 6 ) with two positive initial masses Mi > 

0, M 2 > 0 (for a.e. 0 < t < T) and 

ii) to dissipates the entropy (11) (and thus equally the relative entropy (13) j 
according to the entropy dissipation functional ( 12 ), i.e. satisfies the fol¬ 
lowing entropy dissipation law 

E{u,v,w)(ti)-\- f D{u,v,w){s) < E(u,v,w){tif), for a.e. 0 < to < ti < T, 
Jto 

(14) 

which implies that the entropy E(u,v,w) of the solution remains bounded 
and thus {u{t),v{t),w{t)) G ((LlogL)(fl))^ for a.e. t G [0,T). 

Then, as long as such a solution (u,v,w) exists (i.e. for all t < T < -\-oo) 
and satisfies i) and ii), it also satisfies the following exponential decay toward 
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equilibrium: 


\\u{t) - MoollLi(n) + lk(i) - '!^oo||Li(n) + ||w(t) - ■WoollLi(n) 

<C{E{uo,vq,wo) - E{uoo,Vao,Woo))e~^*, a.e.0<t<T, (15) 

for a constant C and a rate K, which can both be estimated explicitly in terms 
of the parameters a, > 1, di,d 2 ,d^ > 0, Mi, M 2 > 0 and the Poincare- and 
the Logarithmic Sobolev constants P(Ll) and L(Ll) of the domain 

Remark 1.2. The weak entropy dissipation law (14) certainly holds as an equal¬ 
ity for all 0 < to < ti < T for any smooth solutions of suitably truncated 
approximations of Systems (7) - (9). Thus, the weak entropy dissipation law 
(14) also holds as an equality for all classical and weak solutions, which satisfy 
(ypv^ — w'^) G pP for a p > 1, since the terms in the entropy dissipation are 
then uniformly integrable, which allows to pass to the limit. The existence of 
such weak solutions in with p > 1 is a consequence of Proposition 2.1, for 
instance, for all N and a + /3 ,7 < 2 (see (19) and it’s discussion in the Proof of 
Theorem 1.2) or for all N and general a,/ 3,7 provided sufficiently close diffu¬ 
sion coeff dents di,d 2 ,d 3 , which satisfy (18) for some p' < 2. For less regular 
solutions of (7) - (9), like weak -solutions and renormalised solutions (see 
e.g. [6, 9]), assuming they allow to pass to the (a.e. pointwise) limit to obtain 
the weak entropy dissipation inequality (14) (e.g. by lower-semicontinuity, ...), 
Theorem 1.1 still yields exponential convergence to equilibrium in (15) (even 
if passing to the limit in a suitable truncated version of the entropy dissipation 

( 12 ) might not be clear). 

Remark 1.3. Note that in cases, where e.g. classical solutions to system (7) 
- (9) can be established, exponential convergence towards equilibrium in higher 
Lebesgue or Sobolev norms can be proven by an interpolation argument between 
the exponential decay of Theorem 1.1 and polynomially growing bounds 
(compare [3]). Depending on the nonlinearities and the space dimension, Sobolev 
norms of any order may be created even if they do not initially exist thanks to 
the smoothing properties of the heat kernel. 

The proof of Theorem 1.1 is based on the so called entropy method, which 
aims to quantify the entropy dissipation ( 12 ) in terms of the relative entropy 

(13) . More precisely, we shall prove a functional inequality, a so called entropy 
entropy-dissipation estimate (see Proposition 3.1 below) of the form 

D{u,v,w) > K{E{u,v,w) - E{uac,Voo,Woc)), 

where K is an explicitly computable positive constant and (u, v, w) are non¬ 
negative functions, which satisfy the conservation laws (6). The convergence 
to equilibrium as stated in the proof of Theorem 1.1 follows then from applying 
this entropy entropy-dissipation estimate point-wise in time as long as solutions 
dissipate entropy according to (14) (and obey the conservation laws ( 6 )). 

Previous results obtained by applying the entropy method to prove explicit 
exponential convergence to equilibrium for reaction diffusion systems were ob¬ 
tained in e.g [3, 4, 7] for systems featuring quadratic nonlinear reactions. In this 
paper, we present new ideas in proving an entropy entropy-dissipation estimate 
for arbitrary (super)-hnear, monomial reactions rates, which includes general 
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nonlinear mass action law models for any reaction all + /3V O 7 W, see e.g. 

[19]. 

As second main result of the paper, we shall prove the following global exis¬ 
tence Theorem, which partially extends the previously known existence results 
of the System (7) - (9), (4) and (5) by applying duality estimates, which were 
recently developed in [ 2 ]: 

Theorem 1.2. Let LI be a sufficiently smooth bounded (dLl G a > OJ and 

connected open set 0 /IR^ (N > 1). Let di,d 2 ,d 3 > 0 be three strictly positive 
diffusivity constants and a,/3,7 > 1 and (uq, vg, wq) G {L°°(Ll))^. Assume that 
\di — dsl (or equivalently \d 2 — dajj is sufficiently small (to be specified below). 

Then, the System (7) - (9), (4) and (5) has a global classical solution. 

Remark 1.4. The result of Theorem 1.2 seems to be new, for instance, in cases 
where a + (3 > ^ are large and the diffusion coefficients are sufficiently close. 

However, up to our knowledge, the problem of global existence still holds 
open cases when the diffusion coefficients are too different from each other and 
2 < ^ < a + (3. 

We have organised this paper in the following manner. In Section 2, we first 
recall previously proven existence results and prove then Theorem 1.2, i.e. the 
global existence of classical solutions for the System (7) - (9) subject to non¬ 
negative bounded initial data (ug, vg,wo) and homogeneous Neumann boundary 
conditions (4) under the assumption that di (or ^ 2 ) is sufficiently close to dg. 

In Section 3, we establish an entropy entropy-dissipation estimate (Propo¬ 
sition 3.1). Finally in Section 4, we prove a Cziszar-Kullback type inequality 
for the relative entropy functional (13) and show, via a Gronwall argument, 
the exponential decay in L^ towards the equilibrium (uoo,Voo,Woo) as state in 
Theorem 1.1. 


2 Existence Theory 

In the following, we denote LIt = (0, T) x Lt and for p G [1, -foo) 

\\u{t)\\Lp{n) = (^J^\u{t,x)f dx'j , \\u\\LP(nT) = dtdx^ , 

ll'u(t)||L~(n) = esssup |u(t,a;)|, = esssup \u{t,x)\. 

For the sake of completeness and for the reader’s convenience, let us recall the 
main known results on global existence of solutions to System (7) - (9) subject 
to the homogeneous Neumann boundary conditions (4) and non-negative initial 
data (5). At first, we shall provide precise definitions of the notions of solutions 
see e.g. [11, 15]: 

By a classical solution to System (7) - (9), (4) and (5) on LIt = (0,T) x Ll, 
we mean a triple of functions {u,v,w) such that (at least) 

i) iu,v,w) G C{[0,T); L\n)^) n L^i[0,T] X G {0,T), 
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ii) Vp G [l,+oo) we have u,v,w G LP{V,t) and Vfc,£ = 1.. .N 

dtu,dtv,dtw G LP{nT), d^^u,dxkV,dxkW G 
^XkXi'^ j G IjP 

iii) the triple (m, v , w ) satisfy the equations (7) - (9) and the boundary condi¬ 
tions (4) a.e. (almost everywhere) on fly and on (0,T) x i90, respectively 
in the sense of traces. 


By a weak solution to System (7) - (9), (4) and (5) on fir, we denote solutions 
essentially in the sense of distributions or, equivalently here, solutions in the 
sense of the variation of constants formula with the corresponding semigroups. 
More precisely, we assume that — u°‘v^ G L^(Ot) and 


u{t) 

v{t) 

w{t) 


Sd^{t)uo + a f Sdi{t 
Jo 


Sd2{t)vo + fJ 



s){w'^{s) — u°'(s)v^(s)) ds 
s){w'^{s) — u“(s)t>^(s)) ds 


Sd 3 {t)uo+j / Sdsit - s){-w^{s)+ u°‘{s)v^{s))ds 


Jo 

where Sd^i-) is the semigroup generated in L^(n) by —diA with homogeneous 
Neumann boundary condition, 1 < z < 3. 


Provided non-negative initial data uo,vq,wo G the local-in-time ex¬ 

istence and uniqueness of non-negative and uniformly bounded solution to (7) 
- (9) are known (see e.g. [16]). More precisely, there exists T > 0 and a unique 
classical solution {u,v,w) of (7) - (9), (4) and (5) on Qt- If Tmax denotes the 
maximal time of existence, then the solution triple {u, v, w) ceases to remain 
bounded in the sense that 


(Tmax < +Oo) ^ lim (llM(t)llLoo(o) ||z;(t) ||loo(o) -F ||w(t)||L=o(0)) = -Foo. 

max 

Moreover, it is well known that in order to prove global-in-time existence (i.e. 
Tniax = +oo), it is Sufficient to obtain an a-priori estimate of the form 

Vt G [0,Tmax), ||u(t)||L=o(n) -f ||u(<)||^oo(n) -f ||■^«(^)||L~(n) < (16) 

where : [0, -|-oo) —>■ [0, -l-oo) is a non-decreasing and continuous function. 

However, an estimate like (16) is far of being obvious for our System except 
when the diffusion coefficients di, ^ 2 , da are equal, i.e di = da = da = d. In this 
very special case, we have indeed that Z = fj'yu + ajv + 2afiw satisfies 

( Zt — dAZ = 0, (0, Tniax) X 11, 

{E)lnix)-V^Z = 0, (0,T„,ax) X an, 

\^z{o,x) = Zo{x), X G n, 

where Zq{x) = l3juo{x)+a'^vo{x)+2aPwo{x). We may then deduce by the max¬ 
imum principle for parabolic operators (which doesn’t hold for general parabolic 
systems) the following L°° a-priori estimate 

Vt < Tmax, l|2'(t)||L~(n) < ||■^o||L“(n)■ 
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Together with the non-negativity of the solution triple {u,v,w), this implies 

ll'u(i)llL“(n) + ||w(t)||L=o(n) + l|ii'(0llL“(n) < ll■^o||L“(n) 

< \\/3jUo + ajVo + 2a/3wo\\L°-{Q), Vt < Tmax, 

and u{t), v{t) and w{t) are uniformly bounded in L°°{V,) for all f > 0 and 
therefore T^ax = +oo. 

In general, the diffusion coefficients in System (7)-(9) are different from 
each other and the lack of a maximum principle makes the question of global 
existence considerably more complicated. 

Previously, the following cases have been studied by several authors: 

1. Case a = /I = 7 = 1: In this case, global existence of classical solutions was 
first obtained by Rothe [16] for dimensions N < 5. Later, global classical 
solutions have first been proved by Pierre [14] for all dimensions N and 
then by Morgan [13], Martin and Pierre [12] and Feng [ 8 ]. Exponential 
decay towards equilibrium has been shown by Desvillettes-Fellner [3] for 
bounded solutions. Global existence of weak solutions has been proved by 
Laamri [10] for initial data mq, vq and wq only in L^{Q). 

2. Second case 7=1 and arbitrary a, [3 > 1: In this case, global existence of 
classical solutions has been obtained by Feng [ 8 ] in all dimensions N and 
more general boundary conditions. It is also included in the results of 
[13]. 

3. Third case a -f /? < 2 or 7 < 2: In this case, Pierre [15] has proved global 
existence of weak solutions for initial data uq, vq and wq in LF’iJl). 

4. Fourth cases a-|-/ 3 < 7 or(l< 7 < and for arbitrary a,/3 > 1). 
Laamri proved in [11] global existence of classical solutions to System 
(7)-(9) , (4) and (5) in the following cases: 

(a) a-f/3 < 7 , 

(b) (di = da or d 2 = da) and for arbitrary a,/ 3,7 > 1 , 

(c) di = d2 and for any (a, / 3 ,7) such that a -I- /3 7, 

(d) 1 < 7 < and for arbitrary a, /3 > 1. 

Theorem 1.2 below proves the existence of global solutions under a different 
set of assumptions and thus covers cases, which were open so far. The proof 
is based on various duality estimates, in particular on a recent improvement of 
certain duality estimates (see [ 2 ]), which allow to obtain a-priori estimates 
(p > 2 ) in cases where one of the diffusion coefficients di or d 2 is sufficiently 
close to da. More precisely, we shall recall the following Proposition from [2]: 

Proposition 2.1 (See Proposition 1.1 in [2]). 

Let Lt be a bounded domain of with smooth (e.g. (7^“'"“, a > 0) boundary 
dfl, T > 0, and p G (2,-|-oo). We consider a coefficient function M := M{t,x) 
satisfying 

Q < a < M{t, x) < b < +00 for {t, x) € fir, 
for some 0 < a < b < -|-oo, and an initial datum uq G LP{n). 


Then, any weak solution u of the parabolic system 


dtu — Ax{Mu) = 0 

on 

Or, 

m(0, x) = uo(x) 

for 

X G ft, 

n(x) ■ VxU = 0 

on 

[0,r] X dft, 


satisfies the estimate (where p' < 2 denotes the Holder conjugate exponent of p) 


I|i‘I|lp(Ot) ^ + b Da,b,p’)T^^^ \\uo\\LPin), 

and where for any a,b > 0, q € (1, 2) 

C .+6 „ 


(17) 


D. 


a,b,q •— 


1-Ca4 


b—a 

■,q 2 


provided that the following condition holds 

b—a 


C. 


-,p' 


< 1 . 


(18) 


Here, the constant Cm,q > 0 is defined for m > 0, q G {1, 2) as the best (that is, 
smallest) constant in the parabolic regularity estimate 

||^x^^||L9(nj,) < Cm,q ||/||L9(nT)i 

where : [0, T] x H —>■ IR is the solution of the backward heat equation with 
homogeneous Neumann boundary conditions: 


dtv -b TO AxV = f 

on 

fir, 

v{T, x) = 0 

for 

X G ft, 

n(x) ■ VxV = 0 

on 

[0,r] X dft. 


We recall that one has Cm,q < oo for m > 0, q € (1,2] and in particular 
Cm ,2 < Note that the constant Cm,q may depend (besides on m and q) also 
on the domain kl and the space dimension N, but does not depend on the time 
T. 


Proposition 2.1 is the basis for the following proof of the existence Theorem 


1 . 2 : 


Proof of Theorem 1.2. 

The proof performs a bootstrap based on a-priori estimates derived on various 
duality estimates [15, 2]. As a preliminary step, we shall assume without loss 
of generality that |di — daj < |c ?2 — dal- Otherwise, we rename u ^ v, di ^ d 2 
and a —>■ /? and vice versa. 

First, we add the first equation (7) and the third equation (9) to obtain 

dti'ju + aw) — Aj;{M(t,x){'ju + aw)) = 0, where M{t,x) ^-, 

'yu -I- aw 

and observe that 


0 < a := min{di,d 3 } < M{t,x) < b := max{di,d 3 } < -|-oo, V(t,x) G Ot- 
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Thus, by applying the improved duality estimates derived in [2] and the non¬ 
negativity of the solutions, we have that u,w € where po can be 

arbitrarily large provided that b — a = |c?i — dal is assumed sufficiently small 
depending on the space dimension N, see (18) above. In particular, since Cm ,2 < 
^ we estimate for (18) that 


b^ 

O g-f-b 9 ^ ^ i 

2 - a + b 


(19) 


for all 0 < a < 6 < 00 . Moreover, an interpolation argument between Cm ,2 and 
Cm, 3/2 shows that Cm,q for q € [3/2, 2] is bounded by a continuous function in 
q and, as a consequence, that there exists always an exponent po > 2 for any 
values of di and da, see [2]. Finally, in all space dimensions, by assuming that 
|di — dal is sufficiently small (depending on the space dimension N), we can find 
an exponent 

Po =Po(|di - dal) > min{ 7 ,a-f ^}, 

such that at least one of the reaction terms on the right-hand sides of (7)-(9) 
is well-defined in L^. 

Next, we observe that 


/3(dtu - diA^u) = a(dtv - d 2 Axv) 

and by applying another duality argument [15, Lemma 3.4], we obtain 
C'ilkllLp(aT) — ll'“llLP(nT) + 1 ^ C'2 (||t'||LP(nT) + l) > & ( 1 , +00). 

Thus, we have that 

u,v,w € LP°{nT), Po =Po(|di - dal) > min{7,a-|-/3}. (20) 

In the following, we aim at bootstraping (20). We begin by observing that 

and G L'^+o (Ot), 

where the later follows from Young’s inequality, i.e. 

„ a -f 

u < ---. 

a + P 

Then, parabolic regularity applied to the equations (7), (8) implies that (see 
e.g. [5]) 

u,v& (Ot), - > ^ (21) 

qi Po N + 2 

Similarly, equation (9) yields 

weL^H^r), (22) 

qi Po N + 2 

By recalling that 

"/{dtu - diA^u) = a{-dtw + dsA^w), "/{dtv - d 2 Aa,v) = l3{-dtw -I- d^A^w) 
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and by applying [15, Lemma 3.4], we have due to the first equality 


< ||'«||LP(nT) + 1 < 1^2 (llti^llLP(nT) + l) I Vp G (1, +oo), 

and analogous 

C’i||w'I|lp(Ot) — ll'*^llLp(nT) + 1 ^ 1^2 (II'u^IIlp(Ot) + l) I Vp € (l,+oo). 

We are thus able to conclude from (21) and (22) and the above norm equiva¬ 
lencies that 


^ \ . mm{7,a-F/3} 2 

u,v,w € LP^inT), —> - .. I o - 

Pi Po iV -I- 2 

and we calculate that pi > po if and only if 

iV-f 2 


(23) 


Po > 


(mm{ 7 , a + P} — 1). 


In this case, iterating the arguments between (20) and (23) constructs a sequence 
Pn+i > Pm which leads after finitely many steps to an index p„ such that 


min{ 7 , a + P} _ 2 ^ ^ 

Pri fV -1- 2 

and a final bootstrap step yields Pn+i = oo. Thus, the solutions u,v,w are 
bounded in L°°(nT) for any T > 0 and the existence of global classical solutions 
follows by standard arguments. □ 


Remark 2 . 1 . We also remark an alternative approach (with a somewhat dif¬ 
ferent condition) to start the bootstrap in Theorem 1.2, which is based on the 
duality estimates of e.g. [15] and was kindly pointed out to us by a reviewer: 
By adding the first eguation (7) and the third equation (9), we obtain 

dtiiu aw) — diAxiju aw) = 0(^3 — di)AxW. 

Then, by the duality technique (see e.g. the proof of [15, Lemma S.)]) it follows 
directly for any p G (1, -|-oo) and some C = C{p, T) that 

||7'u-haw||LP(nT) < C* [l + Mi - d3||k||LP(nT)] • 

Thus, with a\\w\\LP< || 7 u-|-au;||iP(nr); the required estimate for \\w\\lp(q,.j.) 
holds provided that \di — d^l is small enough. 

3 Entropy entropy-dissipation estimate 

In this section, we prove Proposition 3.1, which details an entropy entropy- 
dissipation estimate for E(u,v,w), D{u,v,w) defined in (11) and (12). The 
proof uses the following technical (but essentially elementary) Lemmata 3.1 
and 3.2: 

At first. Lemma 3.1 establishes a kind of entropy entropy-dissipation esti¬ 
mate in the special case of spatially-homogeneous non-negative concentrations 
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satisfying the conservation laws (6). One could also say that Lemma 3.1 proves, 
what would be the key step of proving an entropy entropy-dissipation estimate 
for the homogeneous ODE system associated to (7)-(9). 

Secondly, Lemma 3.2 generalises the estimate of Lemma 3.1 to spatially- 
inhomogeneous concentrations. 

Notations: For notational convenience, we introduce capital letters as a short 
notation for square roots of lower case concentrations and overline for spatial 
averaging (remember that |0| = 1), i.e. 

U — , U QQ — ; U — / U (t, x) dx, 

Ja 

and analog for V and W. Finally, we denote by ||E ||2 = \F{x)\^ dx the square 

of the norm for a given function E : D —>■ R. 


We begin with; 

Lemma 3.1. Let Uoo, Voo, Woo denote the positive square roots of the steady 
state Uoo,Voc,Woc- Let a,b,c be non-negative constants satisfying the conserva¬ 
tion laws (6), i.e. 

-I- ac^ = Mi = -|- aW^ and 76 ^ -|- = M 2 = -j- fiW^. 

Then, 

(a - Uoo? + {b- Voo? + (c - Woo? < C (a“ b^ - c<? , (24) 

where C = ( 7 ( 0 , /3, 7 ,Mi,M 2 ) is an explicitly computable positive constant. 

Proof of Lemma 3.1. The proof exploits a change of variables introducing 
perturbations pi,fi 2 ,P 3 around the equilibrium values t/ 00 ,14o, Woo, i-e. 

a = Uoo? + hi), b = Voo{l + h 2 ), c = Woo{l + hs)- ( 25 ) 

Since the constants a, b, c are non-negative, the new variable are bounded below, 
i.e. hi,h2,h3 > -1- 

Moreover, the conservation laws (6) rewrite in the new variables into the 
relations 


jUlhi{2 + pi) = -aW^h3{2 + p3), (26) 

jV?h2i2 + p2) = -PWlh3i2 + h3). (27) 

Since the mapping p 1 —>■ p{2 -\- p) is strictly monotone increasing on [—l,oo), it 
is straightforward to show that solving (26) and (27) in terms of p^ defines two 
strictly monotone decreasing functions p^ 1 —>■ piips) and p^ 1 —>■ p 2 {p 3 ), which 
have a unique zero crossing pi{0) = 0 = P2{0) at p 3 = 0. More precisely, 

hiihs) = \J^- 7^^3(2 -h P 3 ) - 1 , Ail = 0 ^3 = 0 , Ai3e[-l,oo), 

h 2 {h‘i) = \J^~ ^^Lz {2 + P 3 ) - 1 , p 2 ?t)^p 3 = Q, P3€[-1,oo). 
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Moreover, the conservation laws ( 6 ) or equivalently the relations (26) and 
(27) imply that the range of admissible /rs is also bounded by above (since 
7*1 (Ms), /^ 2 (M 3 ) are real and bounded below): 


Given these properties, we are able to express and 1 x 2 in terms of ^3 in the 
following way: 

Ml = RlifJ-s) '■= ^jjT 

^i 2 = -R2{^l3)^^3, R2{^l3) ■= 

and it is straightforward to check that the functions i?i(/i 3 ) and i? 2 (M 3 ) are 
strictly monotone increasing and satisfy the following positive lower and upper 
bounds: 

0 < i?l(-l) < = i?l( 0 ) < Rl{^J. 3 , 7 nax) < 00, 

0 < i?2(—1) < = ^^ 2 ( 0 ) < R2{^J'3,max) < OO. 


Next, we insert the ansatz (25) into (24), use the identity and 

observe that in order to prove Lemma 3.1, we have to show that 

^ ^ 

Wc 2 ((1 + Ms)"*' ~ (1 + Mi)“(l + ^2)^) 

for a constant C. We remark that for arbitrary perturbation /ri,/r 2 ,/r 3 > — 1 
the fraction on the left-hand side of (28) will not be bounded as there are many 
more zeros of the denominator than for the nominator, which vanishes only at 
Ml = M 2 = M3 = 0 . 

However, the conservations laws restrict the admissible perturbations mi = 
-i?i(M 3 )M 3 and M 2 = -MJ 2 (M 3 )M 3 for M 3 G [-1, M 3 .maa:]- Thus, we estimate the 
numerator of (28) 


u'Ly:{ + vliil + < fiX 


TJ 2 T /2 

l + ^Rl+ 


'^*00 


^ 00 


R 2 


<C^ll (29) 


for a constant C = C{a, P, 7 , Mi, M 2 ). 

The denominator of (28) can be estimated in the following way: We assume 
hrst that ^3 < 0, which is equivalent to ^i(/i 3 ) > 0,M2(M3) > 0 by the mono¬ 
tonicity properties of and M 2 (M 3 )- Thus, we have (1 -I- ms)^ < (1 + Ms) 

and (1 -I- /ii)“ > (1 -I- Ml) and (1 -I- M2)^ > (1 + M2)- Altogether in this case, we 
have 


|(1 + M 3 )^ - (1 +Mi)“(l + M2)^| > (1 + Mi)(l + M2) - (1 +M3) 

= Ml + M2 + M1M2 - M 3 > Ml - M 3 = (Ri + 1) ImsI (30) 

since mi > ~1 implies M2 + M1M2 > 0. 
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In the case Hs > 0 and thus ^ 1 (^ 3 ) < 0 ,^ 2 (M 3 ) < 0, we have (1 + ^ 3 )'^ > 
(1 + ^ 3 ) and (1 + /ii)“ < (1 + ^ 1 ) and (1 + ^ 2 )^ < (1 + M 2 ) and estimate in a 
similar way 

|(1 + ^ 3 )^ - (1 + Mi)“(l + M 2 )^| > (1 + Ms) - (1 + Mi)(l + M 2 ) 

= ^3 - /ri - ^2 - M 1 M 2 > Ms - Ml = (1 + ^ 1 ) ImsI (31) 

since Mi ^ implies —^2 — M 1 M 2 > 0- 

Altogether, by (30) and (31), we estimate the denominator of (28) by 

((1 + Ms)^ - (1 + Mi)“(l + M2)^)' > fil , 

which proves with (29) the statement of the Lemma 3.1. □ 

The following lemma extends Lemma 3.1 to non-negative functions U, V, W 
whose squares C/^, V^,W^ satisfy the conservation laws (6): 

Lemma 3.2. Let (Uca, Voo, IToo) denotes the positive square roots of the steady 
state {uacyVooiWoo) o,nd [/, V, IT : 17 —>■ IR fee measurable, non-negative functions 
satisfying the conservation laws (6), i.e. 

-flP + aW = Ml = -h aW^, jW -h = M2 = 

Then, the following estimate holds for any a,P,'f > 1 

\\U - UooWl + 111/ - Tooll^ + IIW^ - l^ooll^ < KiWW^ - U’^v^wl 

+ K2 (lie/ - U\\l + 11^ - ^Il2 + - ^ 112 ) , (32) 

for various constants Ki and K2 depending only on a,P,j > 1 , di,d2,d3 > 0 
and Ml, M2 > 0 . 

Proof of Lemma 3.2. 

Step 1 : In order to prove (32), we shall first show that 

||^7_t/«y/3||2 > 1 ^ 

- C {\\U -U\\l + \\V -V\\l + \\W -W\\l) , (33) 

for some constants C. We remark that due to Jensen’s inequality and the 
conservations laws (6), we have the following natural bounds for the averages 

U<VW<Cm, V<VV^<Cm, W <VW <Cm, 

for a constants Cm = CM(a,/3, 7 , Ali, M 2 ). Thus, also the reaction rate 

\U"V^ -W^\< Cm, 

is bounded in terms of the (positive) initial masses Mi, M 2 by various constants 

Cm = Cm{,oc, P, 7, Ml, M2). 

Moreover, we introduce the following deviations with zero mean value: 

6 i{x) = U-U, 62 ix) = V-V, 63{x) = W-W, Sf = S^ =6^ = 0. (34) 
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Step la: We consider for a constant K > 0 the set 


5:={xer! I |(5i| <i^, 1^21 \S 3 \<K}. 

The set S constitutes the part of the state space (C/, V, W) = {U + 6 i,V + 
S 2 , W + 153) which includes the states close to equilibrium at which ([/, V, W) = 
{Uoo^Voo.Woo) and ((5i,(52,<53) = 0 holds. Using the bounds —U < 61 < K, 
— V<S 2 <K and —W<53<K (where to lower bounds are always satisfied 
due to definition (34)), we have with a,/3,7 > 1 

= {U + 5i)“(F + S 2 f = + Ti{x){ 6 i + S 2 ), 

WT- = (W + Ssy =W + T2 {x)S3, 

for bounded remainder terms Ti and T2, which satisfy |Ti(x)| < 

and |T2(a:)| < C{K'^~^), respectively, for all x G S. Thus, we can factorise and 

estimate with Young’s inequality 

_ wy\y^s) = \\u"v^-w^Wh^s) 

+ 2 j {U°V^ -W^)(TySi + S 2 ) - T253) dx 
+ ||T'i((5i + S 2 ) - T2(53||^2(s) 

> 11^2(5) - ||Ti(,5i +^2) -T 253 11^2(5) 

> 1 ||c7“F^ - wy\y^s) - CiK, Ml, M2)(^ + 1 + ^), 

(35) 


where we have used that |Ti| < C{K), IT 2 I < C{K). 

Step lb: It remains to consider the orthogonal set of all points 


5'-^ := {x e fl I |(5i| > K, or |52| > K, or 1^31 > K} . 
By Chebyshev’s inequality it follows that 






\{Sl>Ky\<^. 


Moreover, we observe that |{|(5i| > K}\ = |{5i > K‘^}\. This implies that 
< C{K){ 6 f + 62 + i5|). Therefore, since < Cm is bounded, 

we have 

||C7V-W^||i2(sr) <C(if,Mi,M2)(^ + l + ^) 

and 


> 0 

> ||I7V-W'^||i2(sr)-C(X,Mi,M2)(^ + ^ + ^). (36) 

Altogether, the estimates (35) and (36) finish the proof of the estimate (33) and 
conclude Step 1. 


15 




We remark that estimate (36) could be restricted to the subset 


S^:= {x€S^\ \U‘^V^-W^\<fx}, 


for any fj, > 0 since the estimate (33) holds naturally at all points x € where 
|Vf^7 _ ^ holds, i.e. 


>min 




\u"v-w'' 




w 


L^(s^)- 


Step 2 : Inserting estimate (33) into (32) allows now to apply Lemma 3.1 after 
expanding the left hand side of (32) in terms of 


Si{x) = U-U, 62 {x) = V-V, 63 {x)=W-W, 

Moreover, we shall adapt the ansatz used in Lemma 3.1 in the following fashion 

[p = c/^(l + Ml)^ v^ = vl{l + ^l2)^ m = wl{l + ^l3f, 

and —1 < which recovers the relations (26) and (27) as well as the 

subsequent study of ^i(/i 3 ) = —Ri{n 3 )n 3 and ^jl 2 {p^ 3 ) = —R 2 i^J■ 3 )^J ■3 as functions 
of /i 3 € [-1, fi 3 ,max] (sce Lemma 3.1). 

The definitions of (34) imply readily for the right-hand side of (32) that 

\\U-U\\1=TP-U" =¥, = \\S,\\l (37) 

and analog and . Thus, — = VlP — U and 

^ yP-sF 

it follows that 


U=U^{l + fii)- ^ _ Sl (38) 

Vu'^ + u 

and analog expressions hold for V and W. Therefore, the first terms on the 
left-hand side of (32) expands with (38) to 

= + (39) 

Vu^ + u 

and analog for \\V — Voo\W and ||tT — lToo||i. Note that (39) is quadratic in the 
perturbation and the deviation (5i. 

We point out that in (39) the expansions in terms of Sf is unbounded for 
vanishing U"^ > . Although such states with ~ 0 or 14^ ^ g or W'^ ~ 0 (in 

which the positive initial mass is concentrated in the remaining non-vanishing 
concentrations) are highly degenerate states far from equilibrium, we are nev¬ 
ertheless forced to distinguish two cases. 

Firstly we consider the 

Case > e^, > e^, > £^\ for a constant e > 0 to be chosen later in the 
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opposite cases for small V^, W^. In this case, the expansion (39) is not 

degenerated and we have , ^ ^ ^ < C{e) and 

y/m+U VW^+W y/W^+W 


2Uoo 2Woo 2W^ 

y/W + u’ y/W + W' y/W + W 


< C(e, a,/ 3 , 7 ,Mi,M 2 ). 


(40) 


Moreover, we estimate first by using (38) that 


||C/V^ [(1 + Aii)“(l + fi2f - (1 + M2)1' 

-C(e,a,/3,7,Mi,M2)(^ + l + ^). (41) 

Next, we insert (39) into the left-hand side of (32) (and use (40)) and insert 
(33), (37) and (41) into the right-hand side of (32). Thus, in order to finish the 
proof of the Lemma, it is sufficient to show that 


uIm! + + Wlfil < K,W^ [(1 + Mi)“(l + - (1 + M2)1' 

+ (/L 2 - c(e, a, /3, 7 , K, iLi, Mi, M 2 )) (^ +1 +1) , (42) 

Now, provided that we chose K 2 > C{e,a, K, Ki, Mi, M 2 ) sufficiently 
large, the inequality (42) follows directly from Lemma 3.1, in particular from 
the reformulated estimate (28), which is sufficient to prove estimate (24). 

Secondly we treat the 

Case or or < e^: In these degenerate cases, the expansion 

(38) is not useful. However, since the corresponding states (If, V, W) are far 
from equilibrium, the inequality in (32) can be established by direct estimates 
(i.e. without Lemma 3.1). 

First, we observe that the left-hand side of (32) is bounded in terms of the 
conservation laws ( 6 ) 


\\u - C/00II2 +11^ - KX.II2 + - w^oolli < 


/ {u + v + w + 

Uqc H“ '^00 + Woo)dx < CM{Mi,M2,a,P,'y). 

Jn 


(43) 


Next, we consider the case that U < e^, where e is a constant to be specified 
below. Rewriting the conservation law 'yU'^ + aW^ = Mi using U'^ — U =61 
and =^, we estimate with U < that 


Thus, since Mi > 0 

(w^ _ 2 W'^U^V^ 

> (Ml _ 2W'V^e°‘ > C > 0, 

\ a ^ a a J 


17 









for a strictly positive constant C > 0 provided that 6 f + i5| is sufficiently small 
and e is subsequently chosen small enough. 

Therefore, there exist two constants Ki and K 2 such that 


\\u - C/00II2 + 11^ - Kx.||i + ||VF - IToolli < Cm 

<Ki(W'*-U"V^y+K2(^ + M + '^y (44) 


In the opposite case that 6 f + S§ is not sufficiently small (44) still holds for a 
sufficiently large constant K 2 . 

The case that V < can be treat analog to the case U < above. 

Finally, in the case that W < we estimate the conservation laws yC/^ + 
aW^ = Ml and = M 2 like above as 


7 


a 

7 


((5| + , 







which yields 


(Tr-t7V)“>(^-5;-^(5|+.^) 

- 2 U°‘V^e^ > C > 0, 




for a strictly positive constant C > 0 provided that 5l, (5|, are sufficiently 
small and e is subsequently chosen small enough. 

Similar to above, this yields (44) and together with (33) finishes the proof 
of the Lemma 3.2. □ 


Next, we recall the Logarithmic Sobolev inequality, which holds on a bounded 
domain (without confining potential) as a consequence of the inequalities of 
Sobolev and Poincare: 

Lemma 3.3 (Logarithmic Sobolev inequality). Let (p : fi —)■ R where Lt is a 
bounded domain in such that the Poincare (-Wirtinger) and Sobolev inequal¬ 
ities hold: 

~ /n ‘^^^111,2(0) — 11^ ‘^llL2(n)i 

llv^llL<!(n) — C'i(n) 11^ ‘P\\'L^{n) + C 2 {LI) ||:/3|||2(f2) > q ~ 2 ~ W^ 

Then, the logarithmic Sobolev inequality 

/ - ] dx < L{n,N)\\W,,(p\\l2(^a) (45) 

Jn \ ll‘^llL 2 (o)y 

holds for some constant L(Ll,N) > 0 in any space dimension N. 

The proof follows an argument of Strook [18], see [7]. 

We are now in position to state the entropy entropy-dissipation estimate for 
the relative entropy E — Eoc defined in (11) and the entropy dissipation D (12), 
which holds for admissible functions regardless if or if not they are solutions (at 
a given time t) of the System (7) - (9): 
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Proposition 3.1. Let LI be a connected open set o/R^ {N > 1) with sufficiently 
smooth boundary dLl such that Poincare’s inequality and the Logarithmic Sobolev 
inequality hold. 

Let u, V, w be measurable non-negative functions from Ll into R"*" such that 
J^iju + aw) = Ml and + /3w) = M 2 . Then, 

D{u,V,w) > K{E{u,V,w) - E{Uoc,Voo,Wca)). 


for a constant K = K{di,d 2 ,d 3 ,a, ffij, Mi, M 2 , P{Ll), L{n)) depending only 
on the positive diffusion coefficients di,d 2 ,d 3 > 0, the stoichiometric coeffi¬ 
cients a,/3,7 > 1, the positive initial masses Mi, M 2 > 0 and the Poincare and 
Logarithmic-Sobolev constants P{Ll) and L{Ll) of the domain LI. 

Proof of Proposition 3.1. 

We begin by rewriting the relative entropy (13) using the following additivity 
property 

E{t) — Eoo = J ^uln^—^ + uln^—^ + wln^—c?x (46) 

+ ulnf—) - {u- Moo) + Mlnf—) - {v - Moo) (47) 

VMoo7 VMoo7 

+ MJlnf ) - {W- MIoo), 

\Woc / 

where the integral (46) corresponds to the relative entropy of concentrations 
{u,v,w) with respect to their averages (u,v,w) and (47) is the relative entropy 
of the averages {u,v,w) with respect to the equilibrium (moo, Uoo, Woo)- 

The integral (46) can be estimated thanks to the Logarithmic Sobolev in¬ 
equality (45) by 


[ u\n(^)dx<L(Ll) [ = 4L(fl) [ dx, (48) 

Jn u 

and similar for v and w. Thus, (46) is bounded above by the Fisher information 
part of entropy dissipation (12). 

On the other hand, the relative entropy (47) can be estimated in the following 
way: We define the function 


T{x,y) 


X ln{x/y) - (x-y) 

iVx-i/y)^ 


Pix/y,!), 


which extends continuously to [0, -boo) x (0, -boo) —R+ and note that thanks 
to the conservation laws (6) the expressions 


(fiiu/uooff) < Cm, Tiv/voo,i) < Cm, t{w/w^,1) < Cm, 

are bounded by a constant Cm {Mi, M 2 ) in terms of the initial masses. Thus, 
by using as in Step 2 of Lemma 3.2 that y/E = \/lP = t/oo(l + yi), we have 

uln[— ] - (m- Moo) < C'm(v^- \/u^)^ < CmU"^ yj, (49) 

^ ^00 

and similar estimates hold for v and W. 
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Next, we estimate the entropy dissipation (12) below by using the elemen¬ 
tary inequality {u°‘v^ — w'^){hi{u°'v^) — ln(rt;'^)) > 4(?7“y^ — Poincare’s 

inequality and the Logarithmic Sobolev inequality and obtain the estimate 

D{u,v,w) > -W'^Wl 

+ ecp {\\u-u\\l + ||v^-F||2 + ||w"-WII2) 

+ 4(l-0)(/ |V,C/|^da;+ / / |V,lP|"da:), (50) 

^Jn Jn Jn ' 

for a constant Cp = Cp{di, ^ 2 , da, P{^)) with the Poincare constant P{ft) and a 
constant 6 € ( 0 , 1 ) to be chosen such that the last term on the right-hand side of 
(50) controls via the Logarithmic Sobolev inequality (48) the first contribution 
to the relative entropy E — i.e. the integral (46). 

Combining the expressions of the relative entropy (13) and the entropy dis¬ 
sipation (12) with the estimates (49) and (50) and choosing 9 in (50) in order 
to control (46) with the last term of (50), it remains to show that 

Cm {Ulfil + Vlfil + Wlfil) < K, \\U^V^ - 1P^||J 

+ K,{\\U-U\\l + \\V-V\\l + \\W-W\\l), 

which is a consequence of Lemma 3.2 after observing that 

lit/ - u^wl < \\u - u\\l + \\u - u^wl = \\u - u\\l + ulfil 

and analog for 11P — Poo 11 2 a-nd 11 IP — IPoo 111- This ends the proof of Proposition 

3.1. □ 

4 Estimates of convergence towards equilibrium 

In this section, we use the estimates of Section 3 in order to prove Theorem 

1.1. We begin with a Cziszar-Kullback type inequality relating convergence in 
relative entropy to convergence in L^. 

Proposition 4.1. For all (measurable) functions u,v,w : ^ ( 1 R+)^, for which 

f^(ju + aw) = Ml > 0 and f^('yv + /3w) = M 2 > 0 holds, we have the following 
Cziszar-Kullback type inequality for the entropy functional E{u,v,w) defined in 
( 11 ).- 

E{u,V,w) - £’(Uoo,t'oo, Woo) 



for a constant C = C{Mi, M 2 , a, /3, 7 ) > 0 depending on the masses Mi, M 2 > 0 
the stoichiometric coefficients a,f3,"f> 1 . 

Proof of Proposition 4.1. 

For a constant Xoc > 0, we define (the continuous extension of) the non-negative 
relative entropy density function 

0 < q{x,Xao) ■= a;ln( —) - {x - Xoo), x e [ 0 ,-|-oo), 
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and rewrite the relative entropy (13) as (see also (47)) 

E — Eac = / uln( —J(ia:+ / r;ln( —)(ia;+ / wlnf—J c?x 
Jq ^uJ Jq \v/ Jq \wy 

+ q{u,Uoo) + q{v,Voo)) + q{w,Wao)- (51) 

In order to Taylor-expand the last three terms in (51), we observe Hrst that 

q'{xoa,Xoa) = 0 , and q"{x,Xoo) = —■ 

X 

Since the non-negativity of the solutions and the conservation laws imply the 
natural a-prior bounds 

_ Ml _ M2 _ 

0 < M, Moo < -j 0 < M, Moo < -, 0 < W, Woo < mUl 

7 7 

Taylor expansion yields the following estimate with 6 G (m, Moo) 

(^(m, Moo) — 9(^007^00) 4“ q iy^oo 7 ^00) (^ '^00)4“f7 (^,Moo) ^ 7 

^ Moo7 (m - Mqo)^ 

- Ml 2 ■ 

Similar estimates hold for q{v,Voo) and q{w,Woo) and we obtain 
q{u, Moo) + q{v7 Voo) + q{w, Woo) 

> C(a,/3,7,Mi,M2) (|m- MooP + |T- MooP + |W - rCooP) 

for a constant C{a, /3,7, Mi, M2) > 0 depending only on a, /?, 7, Mi and M2. 

Secondly, considering the integral terms on the right-hand side of (51), we 
estimate with the classical Cziszar-Kullback-Pinsker inequality (Cf. [17]) 

^Mln(|)dx>^||M-M||i.(^), 

and analog for v and w, for which we have again m,m,mJ < C'(a,/3,7, Mi, M2). 

Altogether, after using Young’s inequality to estimate ||m — MooUlqQ) < 
2 ||m — m||^ 1 (-q) + 2 |m — MooP, we obtain 

E — Eoo > C ^||m — Moolliqa) + ||n — '^^oolliqn) + ||w — > 

for a constant C = C(q;,/3, 7, Mi, M2). This ends the proof of Proposition 

4.1. □ 


a ’ /? /’ 


We now are in a position to state the proof of Theorem 1.1. 

Proof of Theorem 1.1. 

The entropy dissipation law (14) yields for the relative entropy with respect 
to the equilibrium, i.e. '0(<) := E(u,v,w){t) — £'(moo, Mqo, m;oo), the following 
inequality for a.e. 0 < to < ti < T 

< i’ito) — f E{s), for a.e. 0 < to < ti < T, 

Jto 
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which rewrites with Proposition 3.1, i.e. D > Ktjj{s) into 

V'(^o) > 4’i'ti) + for S'-®- 0 < to < ti < T. (52) 

Jto 

We can now apply a Gronwall argument as stated in [1] (see also [20] for the 
generalisation where inequality (52) holds only almost everywhere for an inte- 
grable function ijj € ([0, T)) and obtain exponential convergence of the relative 

entropy ip, i.e. 

ipiti) < '0(fo) for a.e. 0 < to < ti < T. (53) 

For the convenience of the reader, we shall recall the proof of (53) in the follow¬ 
ing: First, we perform in (52) the change of variables t = —r and ip{—r) = ip(r) 
and obtain 

/ ro 

ip{s), for a.e. —T< —ti < rp < 0. (54) 

-ii 

Then, we define 4'(r) = J^^^tpir) and calculate with ^(r) = ip{r) > ip{t\) + 
iF4'(r) the well defined derivative 

^ (4'(r)e-^('’+*i)) > (V’(ti) + K'i!{r)) 

Then, integration over [—ti,ro] and division by yields 

^(ro) > ^ ^ 

and further with (54) and ip(r) = 4'(ro) 

^(ro) > ip(ti) 

Then, return to the original variables tg = —rg and ip(—ro) = ip(rg) yields (52) 
and thus by setting to = 0 

E(u,V,w) - E(Uoo,Voo,Woo) < lE(uo,Vo,Wg) - E(Uoo, Voo, Woo)] e~^K (55) 

Finally, the estimate (15) follows from (55) by applying the Cziszar-Kullback 
type inequality in Proposition 4.1. □ 
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